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ABSTRACT: 

We analyze several families of two-dimensional quantum random walks. The feasible region (the 
region where probabilities do not decay exponentially with time) grows linearly with time, as is the 
case with one-dimensional QRW. The limiting shape of the feasible region is, however, quite different. 
The limit region turns out to be an algebraic set, which we characterize as the rational image of a 
compact algebraic variety. We also compute the probability profile within the limit region, which is 
essentially a negative power of the Gaussian curvature of the same algebraic variety. Our me thods 



are based on analysis of the space-time generating function, following the methods of PW02| |. 
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1 Introduction 



Quantum random walk, as proposed by [ADZ9a |. describes the evolution in discrete time of a 
single particle on the integer lattice. The Hamiltonian is space- and time- invariant. The allowed 
transitions at each time are a finite set of integer translations. In addition to location, the particle 
possesses an internal state (the chirality), which is necessary to make the evolution of the location 
no ndetermini stic. A rigorous mathematical analysis of this system in one dimension was first given 
by ABNj^Ol|. The particle moves ballistically, meaning that at time n, its distance from the origin 
is likely to be of order n. By contrast, the classical random walk moves diffusively, being localized 
to an interval of size y/n at time n. 

A very similar process may be defined in higher dimensions. In particular, given a subset 
_E C Z'' with cardinality k and a k x k unitary matrix U, there is a corresponding space- and 
time-homogeneous QRW in which allowed transitions are translations by elements of E and evo- 
lution of chirality is governed by U. When E is the set of signed standard basis vectors we 
call this a nearest neighbor QRW; for example in two dimensions, a nearest neighbor walk has 
E — {(0, 1), (0, — 1), (1, 0), (— 1, 0)}; a complete construction of quantum random walk is given in 
Section 12.11 below . Publish ed work on quantum random walk in dimensions two and h igher began 



around 2002 (see |MBSS02l |). Most studies, including the most recent and broad study |WKKK08l |. 
are concerned to a great extent with localization; this phenomenon is not generic in quantum ran- 
dom walk models and not present in the models we discuss below. The analyses we have seen 
range from analytic derivations without complete proofs to numerical studies. As far as we know, 
no rigorous analysis of two-dimensional QRW has been published. The question of describing the 
behavior of two-dimensional QRW was brought to our attention by Cris Moore (personal commu- 
nication). In the present paper, we answer this question by proving theorems about the limiting 
shape of the feasible region (the region where probaiblities do not decay exponentially with time) for 
two-dimensional QRW, and by giving asymptotically valid formulae for the probability amplitudes 
at specific locations within this region. 

Our analyses begin with the space-time generating function. This is a mu ltivaria te rational 
function which may be derived without too much difficulty. The companion paper |BP07l | introduces 
this approach and applies it to an arbitrary one-dimensional QRW with two chiralities (fc = 2) . This 
approach allows one to obtain detailed asymptotics such as an Airy-ty pe limit in a scaling window 
near the endpoints. As su ch, it im proves on the analysis of ABN^Ol] but not on the more recent 
and very n ice ana lysis of In one dimension, when the number of chiralities exceeds two, 

N. Konno [Kon05l | found new behavior that is qualitatively different from the two-chirality QRW. 



Forthcoming work of the last author with T. Greenwood uses the generating function approach 
to greatly extend Konno's findings. The generating function approach, however, pays its greatest 
dividends in dimension two and higher. This approach is based on recent results on asymptotics 
of multivariate rational generating functions that allow nearly automatic transfer from rational 
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generating functions to asymptotic formulae for their coefficients PW02I : |PW04l : iPWOSt |BP08[. 
Based on these results, analyses of any instance of a two-dimensional QRW becomes relatively 
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(a) limit (b) Exact probabilities at time 200 

Figure 1: Fixed-time empirical plot versus theoretical limit 



easy, although in some cases new versions of the results under weaker hypotheses were required. 
Empirically computed probability profiles such as are shown in figure llbl are explained by algebraic 
computations, leading to limit shapes as shown in figure [Tal We computed probability profiles for a 
number of instances of two-dimensional QRW. The pictures, which appear scattered throughout the 
paper, are quite varied. Not only did we find these picture s visually intriguing, but they pointed us 
toward some refinements of the theoretical work in [PW02|, which we now describe, beginning with 
a more detailed description of the two plots. 

On the right is depicted the probability distribution for the location of a particle after 200 steps 
of a quantum random walk on the planar integer lattice; the particular instance of QRW is a nearest 
neighbor walk {E — {(0, 1), (0, —1), (1, 0), (—1, 0)}) whose unitary matrix is discussed in Section^ 
Greater probabilities are shown as darker shades of grey. The feasible region, where probabilities are 
not extremely close to zero, is the diamond with vertices at the midpoints of the 400 x 400 square. 
The feasible region appears to be a slightly rounded diamond. 

In his Masters Thesis, the second author computed an asymptotically valid formula for the 
probability amplitudes associated with some instances of QRW. As n ^ oo, the probabilities become 
ex ponent ially small outside of a certain algebraic set H, but are Q{n~^) inside of H. Theorem 4.5 
of BraOTj l proves such a shape result for a different instance of two-dimensional QRW and conjectures 



it for this one, giving the believed characterization of S as an algebraic set. The plot in figure [Tal 
is a picture of this characterization, constructed by parametrizing S by patches in the flat torus 
To :— (]R/27rZ)^ and then depicting the patches by showing the image of a grid embedded in the 
torus. 

When the plot was constructed, it was intended only to exhibit the overall shape. Nevertheless, 
it is visually obvious that significant internal structure is duplicated as well. Identical dark regions 



2 



in the shape of a Maltese cross appear inside each of the two figures. To explain this, we consider 
the map <i> : T ^ whose image produces the region S. Let V denote the pole variety of the 
generating function F for a given QRW, that is, the complex algebraic hypersurface on which the 
denominator H oi F vanishes. Let Vi denote the intersection of V with the unit torus T. It is easy 
to solve for the third coordinate z as a local function of x and y on Vi and thereby obtain a piecewise 
parametrization 



of Vi by patches in . Theorem 13.31 extends the results of PW02| | to show that each point z of Vi 
produces a polynomially decaying contribution to the probaiblity profile for movement at velocity 
(r, s) which is the image of z under the logarithmic Gauss map n of the surface Vi at z: 

f dH dH dH 
n(z) := [-^,y^,-^ 

Formally, the n maps into the projective space MP^, but we map this to by taking the projection 
7r(r, s,t) :— {r /t, s/t,l). In other words, the plot is the image of the grid (Z/IOOZ)^ under the 
following composition of maps: 

(Z/100Z)2 ^S'xS^ V ^ RP2 ^ _ 

The intensity of an image of a uniform grid of dots is proportional to the inverse of the Jacobian 
of the mapping. The Jacobian of the composition is the product of the Jacobians of the factors, 
the most significant factor being the Gauss map, n. Its Jacobian is just the Gaussian curvature 
(in logarithmic coordinates). The darkest regions therefore correspond to the places where the 
curvature of logVi vanishes. Alignment of this picture with the empi rical am plitudes can only 



mean that the formulae for asymptotics of generating functions given in PW02| | blow up when the 
Gaussian curvature of log Vi vanishes. This observation allowed us to produce new expressions for 
the quantities in the conclusions of theorems in |PW02] , where lengthy polynomials were replaced 
by quantities involving Gaussian curvatures. 

To summarize, the purpose of this paper is twofold: 



In Theorem 14.91 we prove the shape conjecture from Bra07[; further instances of this are 
proved in Theorems 14.21 and 14.71 



In Theorems 13.31 and 13.51 we reformulate the main result in [PW02| | to clarify the relation 
between the asymptotics of a multivariate rational generating function and the curvature of 
the pole variety in logarithmic coordinates. 

The organization of the remainder of this paper is as follows. Section [2] gives some background on 
quantum random walks, notions of Gaussian curvature, amoebas of Laurent polynomials, the mul- 
tivariate Cauchy formula, and certain standard applications of the stationary phase method to the 
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evaluation of oscillating integrals. Section [3] contains general results on rational multivariate asymp- 
totics that will be used in the derivation of t he QRW limit theorems. In particular, Theorem 13.31 
gives a new formulation of the main result of |PW02| | , while Theorem 13.51 proves a version of these 
results in situa tions where the geometry of Vi is more complicated than can be handled by the meth- 
ods of |PW02| |. Finally, Section 2] applies these results to a collection of instances of two-dimensional 



nearest neighbor QRW in which the unitary matrices are elements of one-parameter families named 
S{p),A{p) and B{p), < p < 1. This results in Theorems g^l [12] and US] respectively. The QRW 
in figure [T] has unitary matrix B{l/2), while the following figures show examples of the S{l/2) and 
A(5/9) quantum random walks. 




(a) limit (b) probabilities at time 200 

Figure 2: the 3(1/2) QRW 




(a) limit (b) probabilities at time 200 

Figure 3: the ^(5/9) QRW 
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2 Preliminaries 



2.1 Quantum random walks 



The quantum random walk is a model for the motion of a single quantum particle evolving in Z'' 
under a time and translation invariant Hamiltonian for which the probability profile of a particle 
after one time step, started from a known location, is uniform on the neighbors. Such a process was 

/C^)} C be 



first constructed in 



ADZ93 |. Let d > 1 be the spatial dimension. Let E = {v'-^-', . . . , 
a set of finite cardinality k. Let J7 be a unitary matrix of size k. The set Z'' x iJ indexes the set 
of pure states of the QRW with parameters k, E and U . Let lA^U denote the operator that sends 
(r,v('')) to (r, //v^-')), that is, it leaves the location unchanged but operates on the chirality by U. 
Let fj denote the operator that sends (r,v(^') to (r + v'^-'\ v*^-')), that is, it translates the location 
according to the chirality and does not change the chirality. The product a ■ (Id ® U) is the operator 
we call QRW with parameters fc, E and U . Let us denote this by Q. 



For 1 < i,j < k and r e 



denotes the amplitude at time n for a particle starting at location in chirality i to be in location 
r and chirality j. Let z denote (zi, . . . , Zd+i) and define 



(2.1) 



which denotes the spacetime generating function for n-step transitions from chirality i to chirality 
j and all locations. Let F(z) denote the matrix (i^*-*'-''')i<ij<fc. Let M denote the diagonal matrix 
whose entries are the monomials {z"" : r G E}. When d = 2 we use (x, y, z) for (zi, Z2, z^) and (r, s) 
for r; for a two-dimensional nearest neighbor QRW, therefore, the notation becomes 



and 



F(^'^")(x,y,z)= ^VF(r,s)x'- 

n.r.s 



M = 



/ X 
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An explicit expression for F may be der ived via an elementary enumerative technique known as 
the transfer matrix method 



Sta97 



GJ83l |. F or d = \ a nd a partic ular choice of U (the Hadamard 



matrix) , this rational function is computed in ABN+Olj . In BPOTl Section 3] , the following formula 
is given for the matrix generating function F: 



F(z) = (/ - Zd+iMUy 



(2.2) 
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The (j, j)-entry of the matrix, F^'^'^\ may therefore be written as a rational function G/H where 



H = det{I - Zd+iMU) . 

The following result is easy but crucial. It is valid in any dimension d> 1. Let denote the unit 
torus in C*. 

Proposition 2.1 (torality). The denominator H of the spacetime generating function for a quantum 
random walk has the property that 

izi,...,zd) eTd and H{z) =0 =^ \zd+i\ = I . (2.3) 



Proof: If {zi, . . . ,Zd) e then M is unitary, hence MU is unitary. The zeros of det(/ — z^+iAft/) 
are the reciprocals of eigenvalues of MU, which are therefore complex numbers of unit modulus. □ 

Proposition 2.2. Let H be any polynomial and let V denote the pole variety, namely the set 
{z : H{z) = 0}. Let Vi := V D T^+i. Assume the torality hypothesis (j2.3p . Let p £ Vi be any point 
for which \JH{p) ^ 0. Then Vi is a smooth d-dimensional manifold in a neighborhood of p. 



Proof: We will show that dH/dzd+i{p) ^ 0. It follows by the implicit function theorem that there 
is an analytic function 17 : ^ C such that for z in some neighborhood of p, if (z) = if and only 
if Zd+i — g{zi, . . . , Zd)). Restricting (zi, . . . , Zd) to the unit torus, the torality hypothesis implies 
Zd+i — 1, whence Vi is locally the graph of a smooth function. 

To see that dH/dzd+i{p) ^ 0, first change coordinates to zj — pj ex.p(i9j) and Zd+i — Pd+i exp(i(T). 
Letting H H o exp, the new torality hypothesis is {9i, . . . , 0d) S and H{9i, . . . ,9d,o') — im- 
plies (T G M. We are given ViJ(O) 7^ and are trying to show that dH/da-{0) 7^ 0. 

Consider first the case d ~ 1 and let 9 := 9i. Assume for contradiction that dH /d(j{0,0) = 
^ dH/d9{Q, 0). Let H{9, cr) ~ J^j k>o ^iM^'''^^ be a series expansion for i7 in a neighborhood of 
(0,0). We have 60.0 = 87^ 61.0. Let i be the least positive integer for which the 60. £ 7^ 0; such an 
integer exists (otherwise -ff (0, cr) = 0, contradicting the new torality hypothesis) and is at least 2 by 
the vanishing of dH/d(j{0, 0). Then there is a Puiseux expansion for the curve {H = 0} for which 
c ~ (^1,0^/^0, f)^^^- This follows from EKsd ] although it is quite elementary in this case: as cr, ^ 0, 
the power series without the (1, 0) and {0,£) terms sums to 0(|f?P + \9a\ + |cr|^+^) = o{\9\ + \a\^) (use 
Holder's inequality); in order for H to vanish, one must therefore have bifi9 + b^^ia^ — o{\9\ + |cr|^), 
from which a ^ {bifl9 /b^xY/^ follows. The only way the new torality hypothesis can now be satisfied 
is if ^ = 2 and 6i.o^/&o/ is always positive; but 9 may take either sign, so we have a contradiction. 

Finally, if d > 1, again we must have &o,...,o.£ 7^ in order to avoid i?(0, . . . , 0, cr) = 0. Let r e M'^ 
be any vector not orthogonal to VH{<d) and let G{9, cr) := H{ri9, rd9, cr). Then dG/d9{0, 0) ^ 
= dG/da{0, 0) and the new torality hypothesis holds for G; a contradiction then results from the 
above analysis for the case d= 1. □ 
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A Hadamard matrix is one whose entries are all ±1. There is more than one rank-4 uni- 
tary matrix that is a constant multiple of a Hadamard matrix, but for some reason the "standard 
Hadamard" QRW in two dimensions is the QRW whose unitary matrix is 



1 



-1 



Had 



-1\ 



1 

-1 1 -1 -1 

-1 -1 1 -1 

V -1 -1 -1 1 y 



IKK04: 



WKKK08II as the "Grover walk" because of its relation to the 



This is referred to by Konno 
quantum search algorithm of L. Grover. Shown in figure[4a|is a plot of the probability profile for the 
position of a particle performing a standard Hadamard QRW for 200 time steps. This is the only 
two-dimensional QRW we are aware of for which even a nonrigorous analysis had previously been 
carried out. On the right, in figure l4bl is the analogous plot of the region of non-exponential decay. 



'mm 



(a) limit (b) exact probabilities at time 200 

Figure 4: Moore's Hadamard QRW 



Another 4x4 unitary Hadamard matrix reflects the symmetries of (Z/ (2Z))^ rather than Z/ (4Z): 

This matrix also goes by the name of 5'(l/2) and is a member of the first family of QRW that we will 
analyze. There is no reason to stick with Hadamard matrices. Varying U further produces a number 
of other probability profiles including the families S(j)),A{p) and B{p) analyzed in Section[4l 



C^Had 



/ 1 1 

-1 1 

1 -1 

V-1 -1 



2.2 Differential Geometry 

For a smooth oricntable hypersurface V C M''+^, the Gauss map n sends each point p G V to a 
consistent choice of normal vector. We may identify n{p) with an element of S"^. For a given patch 
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P C V containing p, let n[P] := Uq^pn{q), and denote the area of a patch P in either V or 5*^ as 
A[P]. Then the Gauss-Kronecker curvature of V at p is defined as 

JC:^lun^. (2.4) 

p^P A[P] ^ ' 

When d is odd, the antipodal map on S'^ has determinant —1, whence the particular choice of unit 
normal will influence the sign of /C, which is therefore only well defined up to sign. When d is 
even, we take the numerator to be negative if the map n is orientation reversing and we have a 
well defined signed quantity. Clearly, /C is equal to the Jacobian of the Gauss map at the point p. 
For computational purposes, it is convenient to have a formula for the curvature of the graph of a 
function from M'' to K. 

Proposition 2.3. Suppose that in a neighborhood of the point p, the smooth hypersurface V C W^^^ 
is the graph of a function h mapping the origin to p; that is, in some neighborhood of the origin, 

V = {(x, t) : r = /),(x)}. Let V := VftYO) and Ti. := dct — — - — (0) ) denote respectively the 

gradient and Hessian determinant of h at the origin. Then the curvature of V at the origin is given 
by 

K= ^ 



The square root is taken to be positive and in case d is odd, the curvature is with respect to a unit 
normal in the direction in which the dependent variable increases. 

Proof: Let X : U C R** — > M'*+^ denote the parametrizing map defined by 

X(u) := (ui, . . . , Ud, h{ui, . . . , Ud)) 

on a neighborhood U of the origin. Let tt be the restriction to V of projection of R''+^ onto the first 
d coordinates, so tt inverts X on U. Define a vector 

, / dh dh ^ 

\ dui ' ' ■ ' ' dud ' 

normal to V at X(u) and let N denote the corresponding unit normal N/|N|. Observe that |N| = 
+ |V/i|^, and in particular, that |N(0)| = \/l + |Vp. The Jacobian of tt at the point p is, up to 
sign, the cosine of the angle between the Zd+i axis and the normal to the tangent plane to V at p. 
Thus 

= \^L^ = IZM = . (2.5) 

iN||e,+ii 1-1 v/nw 

The Gaussian curvature at the point p is, by definition, the Jacobian of the map N o tt at p. 
Using J to denote the Jacobian, write N as | • | o N and apply the chain rule to see that 

/C = J(7r(p)) • J(N)(0) • J(| • I )(N(0)) = ^ • J(N)(0) • J(| • | )(V, -1) . (2.6) 
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Here, | • | is considered as a map from M'' x {—1} to S'^: at the point y, its difFcrential is an orthogonal 
projection onto the plane orthogonal to (y, —1) times a rescaling by |(y, — 1)|^^, whence 

^(1 • I )(y) = VTTW'V^TW' ■ (2.7) 

Because N maps into the plane z^^i = — 1 we may compute J(N) from the partial derivatives 
dN,/dxj = d^h/dx^dxj, leading to J(N)(0) = H. Putting this together with (^7)) gives 

^(N)(0) = — (2.8) 



and using (|2.6p and p.5p gives 



rd+2 



proving the proposition. □ 

We pause to record two special cases, the first following immediately from V/i(0) =0. If Q is 
a homogeneous quadratic form, we let \\Q\\ denote the determinant of the Hessian matrix of Q; to 
avoid confusion, we point out that the diagonal elements an of this matrix are twice the coefficient 
of xf in Q. The determinant will be the same when the coefficients of \\Q\ \ may be computed with 
respect to any orthonormal basis. 

Corollary 2.4. Let V be the tangent plane to V at p and let v be a unit normal. Suppose that V is 
the graph of a smooth function h over V , that is, 

V = {p + u + /i(u)v : u e C/ C -p} . 

Let Q be the quadratic part of h, that is, h{u) = Q{u) + 0(|up). Then the curvature of V at p is 
given by 

ic = \\Q\\- 

□ 

Corollary 2.5 (curvature of the zero set of a polynomial). Suppose V is the set {x : -ff(x) = 0} 
and suppose that p is a smooth point of V, that is, \7H{p) ^ 0. Let V and Q denote respectively 
the gradient and quadratic part of H at p. Let Q± denote the restriction of Q to the hyperplane V_l 
orthogonal to V. Then the curvature of V at p is given by 

^ = ^ • (2-9) 

Proof: Replacing H by |V|^^i/ leaves V unchanged and reduces to the case \VH{p)\ = 1; we there- 
fore assume without loss of generality that |V| — I. Letting uj^ + A(u)V denote the decomposition 
of a generic vector u into components in (V) and Vj^, the Taylor expansion of H near p is 

+ u) = V • u + Qj.(u) + R 
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where R = 0(|uj^p + |A(u)||uj^|). Near the origin, we solve for A to obtain a parametrization of V 
by V^: 

A(u) =g^(u) + 0(|u|3). 
The resuh now follows from the previous corollary. □ 



2.3 Amoebae and Cauchy's formula 

Let F = G/H he a quotient of Laurent polynomials, with pole variety V := {z : H{z) ~ 0}. Let 
Log : (C*)''+i denote the log-modulus map, defined by 

Log(z) := (log|zi|, . . . ,log|zd+i|) . 

The amoeba of H is defined to be the image under Log of the variety V. To each component B 
of the complement of this amoeba in M''+^ corresonds a Laurent series expansion of F. When F is 
the {d+ l)-variable spacetime generating function of a d-dimensional QRW, we will be interested in 
the component Bq containing a translate of the negative Zcj+i-axis; this corresponds to the Laurent 
expansion that is an ordinary series in the time variable and a Laurent series in the space variables. 
For QRW, the point is always on the boundary of Bq. In general, all compone nts of the complement 



of any amoeba are convex. For further details and properties of amoebas, see GKZ94I . Chapter 6] 



For any r G M''+^, let f denote the unit vector r/|r|. Two important hypotheses that will be 
satisfied for QRW are as follows. 

The function r • x is maximized over Bq at a specified point ; (2-10) 

we will be primarily concerned with those f for which this maximizing point is the origin, and we 
denote by K the set of r for which this holds: thus for f S K and x G Bq, r • x < with equality 
when X = 0. Secondly, we assume that the set W = W(r) of z = exp(x + iy) such that 

H{2.) = and Viogi?(z) || f (2.11) 

is finite. The set W(r) depends on r only through f. The gradient oi H o exp at the point z e W is 
equal to {zidH/dzi, . . . , Zd+idH/dzd+i) and will be denoted Viogi?(z). It is immediate from (|2.1ip 
that Viogi?(z) is a multiple of the real vector r for all 1 < j < fc. 

Before we proceed we point out a condition under which (|2.1ip is always satisfied. Suppose that 
Vi is smooth off a finite set W, and we let r be some direction such that hypothesis (|2.1ip fails. The 
set W(r) is algebraic, so if it is infinite it contains a curve, which is a curve of constancy for the 
logarithmic Gauss map. This implies that the Jacobian of the logarithmic Gauss map vanishes on 
the curve, which is equivalent to vanishing Gaussian curvature at every point of the curve. Thus, if 
we restrict r to the subset of Vi where K, ^ 0, then hypothesis (|2.1ip is automatically satisfied. 
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The coefficients of the Laurent series corresponding to Bq may be computed via Cauchy's 
integral formula. Define the flat torus To := (M/ (2n'Z))'^^^ . The following proposition is well known. 

Proposition 2.6 (Cauchy's Integral Formula). For any u interior to Bq, 

Or = j exp(-r • u) y" exp(-ir • y)F o exp(u + iy) dy . (2.12) 

Corollary 2.7. Let A := A(r) := sup{f • x : x G Bq}. For any X' < X, the estimate 

|ar'|=o(exp(-A'|r'|)) 
holds uniformly as r' —> oo in some cone with r in its interior. 

Proof: Pick u interior to Bq such that r • u > A'. There is some e > and some cone K with 
r in its interior such that r' • u > A' + e for all r' e K. The function F is bounded on the torus 
exp(u + iy), and the corollary follows from Cauchy's formula. □ 

Note: We allow for the possibility that hypothesis (|2.1ip holds for no points with modulus 1. In 
the asymptotic estimate p.6p below, the sum will be empty and we will be able to conclude that 
flr — 0(|r|^(''+^^/^), as opposed to Odrj"''/^) in the more interesting regime; we will not be able 
to conclude that Or decays exponentially, as it does when r ^ K. This will correspond to the case 
where in fact r e K \ K. Observe also that the finiteness hypothesis (|2.10p is not required for this 
result. 

2.4 Oscillating integrals 

Let M be an oriented d-manifold, let (f> : Ai ^ M. he a smooth function and let ^ be a smooth d-form 
on A4. Say that G Al is a critical point for (p if d(j){p^,) = 0. Equivalently, in coordinates, is 
critical if the gradient vector V(/)(p*) vanishes. At a critical point, (f>(jj) — (t>{p*) is a smooth function 
of p which vanishes to order at least 2 at p = . Say that a critical point for is quadratically 
nondegenerate if the quadratic part is nondegenerate; in coordinates, this means that the Hessian 
matrix 

/ i<i,j<k 



has nonzero determinant. It is well known (e.g., |BH86l : IWon89l | that the integral exp{iX(j){y))A{y) dy 



can be asymp totical ly estimated via a stationary phase analysis. The following formulation is 



adapted from [Ste93|. 



If p 1-^ (xi, . . . , Xd) is a local right-handed coordinatization, we denote by ri[p, dx] the value A{p) 
for the function A such that rj = A(jj) dx. If the real matrix M has nonvanishing real eigenvalues, we 
denote a signature function a{M) := n+{M) — n-{M) where n+{M) (respectively n_(M)) denotes 
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the number of positive (respectively negative) eigenvalues of M . Given (p and 77 as above, and a 
critical point p.^, for 0, we claim that the quantity T defined by 

^(0,ry,p,) :=e-*""/4|detH(0;p,)r^/'??[p,,dx] (2.13) 

does not depend on the choice of coordinatization. To see this, note that the symmetric matrix 
7i has nonzero real eigenvalues, whence ili. has purely imaginary eigenvalues and the quantity 
e~"'^/^| det7i((/);p*)|~^/^ is a —1/2 power of det(i7i), in particular, the product of the reciprocals 
of the principal square roots of the eigenvalues. Up to the sign choice, this is invariant because a 
change of coordinates with Jacobian J at divides ?7[p*,c?x] by J and 7i((/);p^) by . Invariance 
of the sign choice follows from connectedness of the special orthogonal group, implying that any two 
right-handed coordinatizations are locally homotopic and the sign choice, being continuous, must 
be constant. 

Lemma 2.8 (nondegenerate stationary phase integrals). Let cj) he a smooth function on a d-manifold 
M. and let rj be a smooth, compactly supported d-form on M. Assume the following hypotheses. 

(i) The set W of critical points of cj) on the support of rj is finite and non-empty, 
(ii) (j) is quadratically nondegenerate at each p^ G W. 



Ther, 



^exp(^A<^)ry=(^y^' ^ e^^^(^*)^(0,,y,p.) + o(A-('^+i)/2) 



(2.14) 



Remarks. The stationary phase method actually gives an infinite asymptotic development for this 
integral. In our application, the contributions of order X^'^l'^ will not ca ncel, in which case p.l4p 



gives an asymptotic formula for the integral. The remainder term (see Ste93|) is bounded by a 
polynomial in the reciprocals of |V0| and detTi and partial derivatives of </) (to order two) and 77 (to 
order one); it follows that the bound is uniform if (\) and 77 vary smoothly with (i) and (ii) always 
holding. 



Proof: Let {Ma) be a finite cover of M. by open sets containing at most one critical point of (/), 
with each Ma covered by a single chart map and no two containing the same critical point. Let \^aS 
be a partition of unity subordinate to {A/'q}. Write 

I ■= exp{iX(j)) T] 

J M 

as -^a where 

la := / exp{iX(j)) r] ipa ■ 



Wo 
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According to Ste93l . Proposition 4 of VIII. 2.1], when TVq contain s no cr itical point of (j> then is 
rapidly decreasing, i.e, /o,(A) = o(A~^) for every N. According to [Ste93l . Proposition 6 of VIII. 2. 3], 



when A/'q contains a single nondegenerate critical point for then, using the fact that ipaip*) — 1, 

\ / j=i 

where rj = ^(x)dx in the local chart map, {/ij} are the eigenvalues of iH, in this chart map, and the 
principal —1/2 powers are chosen. Summing over a then proves the lemma. □ 

As a corollary, we derive the asymptotics for the Fourier transform of a smooth d-form on an 
oriented d-manifold immersed in K'^+^. Let A4 be such a manifold and let IC{p) denote the curvature 
of M at p. If ry is a smooth, compactly supported d-form on M, denote rj[p] = r][p, dx] with respect 
to the immersion coordinates, and define the Fourier transform t) by 

Jm 

Corollary 2.9. Let K be a compact subset of the unit sphere. Assume that for f e K , the set W 

of critical points for the phase function f • x is finite (possibly empty), and all critical points are 
quadratically nondegenerate. For x € W, let r(x) denote the index of the critical point, that is, 
the difference between the dimensions of the positive and negative tangent suhspaces for the function 
f • X. Then 

f,{r)^C^y'^ ^ e^----r;[x,]^(x,)-i/2e-"^(''-)/4 + 0(A-('^+i)/2) 
^' x.ew 

uniformly as |r| —> oo with r G K . 



Proof: Plugging (/> = r • x into Lemma [2?8l and comparing with p.l3p we see that we need only to 
verify for each x* G W that 

e-"^/^ |detH(0;xOr'/' r][x,, dx] = f^[x,] \IC{x,)\-^/^ ^-^-r(^-^)l^ . 

With the immersed coordinates, a ~ t, and this amounts to verifying that | dct7i((/); x*)| = |/C(x*)|. 
Let V denote the tangent space to M. at x* and let u\, . . . ,Ud be an orthonormal basis for V . Let 
V be the unit vector in direction f , which is orthogonal to V because x* is critical for 0. In this 
coordinate system, express as a graph over V . Thus locally, 

A4 = {x* + u + h{vL)v : u e T'} 

for some smooth function h with /i(0) and V/i(0) vanishing. Let Q denote the quadratic part of h. 
By Corollary [Ml we have /C(x*) = ||Q||. But 

(/)(x* + u + h{vL)v) = (/>(x*) + /i(u) 

whence 7^(0; x*) = Q, completing the verification. □ 
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3 Results on multivariate generating functions 



In this section, we state general results on asymptotics of coeffic ients of rational multivariate gen 



crating functions. These results extend previous work of jPW02l | in two ways: the hypotheses are 
generalized to remove a finiteness condition, and the conclusions are restated in terms of Gaussian 
curvature. Our two theorems concern reductions of the (c?+l)-variable Cauchy integral to something 
more manageable; the second theorem is an extension of the first. 

We give some notation and hypotheses that are assumed throughout this section. Let F — G/H 
be the quotient of Laurent polynomials in d + 1 variables z :— (zi, . . . , Zd+i) and let Bq be a 
component of the complement of the amoeba of H containing a translate of the negative z^+i-axis 
(see Section 1^31) . Assume G OBq and let F — J2r '^rz'^ be the Laurent series corresponding to Bq. 
Let V denote the set {z e C''+^ : H{z) = 0} and Vi := V n T denote the intersection of V with the 
unit torus. Let W :— Vi O {z : WH{z) = 0} denote the singular set of Vi. Let K :— K(0) denote 
the cone of f for which the maximality condition (|2.10p is satisfied with x» = and let Af be any 
compact subcone of the interior of K such that (|2.1ip holds for r E Af (finitely many critical points). 



3.1 When V is smooth on the unit torus 

We start with the definition/construction of the residue form in the case of a generic rational function 
F — P/Q with singular variety Vq . 

Proposition 3.1 (residue form). There is a unique d-form i], holomorphic everywhere VQ does 
not vanish such that rj A dQ ^ Pdz. We call it the residue form for F on Vq and denote it by 
RES (Fdz). 

Remark. To avoid ambiguous notation, we denote the usual residue at a simple pole a of a univariate 
function / by 

residue(/; a) = lim (z — a)f{z) . 

z — ya 

Proof: To prove uniqueness, let r]i and 772 be two solutions. Then (771 — 772) A dQ = 0. The inclusion 
L : Vq induces a map l* that annihilates any form ^ with ^ A dQ = 0. Hence rji = 772 when 

they are viewed as forms on Vq. 

To prove existence, suppose that {dQ /dzd+i){z) 7^ 0. Then the form 

is evidently a solution. One has a similar solution assuming dQ/dzj is nonvanishing for any other 
j. The form is therefore well defined and nonsingular everywhere that VQ is nonzero. □ 
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From the previous proposition, RES (F dz) is holomorphic wherever Vi? ^ 0, and in particular, 
on Vi \ W. 

Lemma 3.2. Let F,G, H,V, Bq,Vi and W be as stated in the beginning of this section. Assume 
torality (|2.3p and suppose that the singular set W is empty. Then may be computed via the 
following holomorphic integral. 

d 



1 

27ri 



r-l 



RES [F dz) 



Vi 



(3.2) 



Proof: As a preliminary step, we observe that the projection tt : V ^ onto the first d coordinates 
induces a fibration of Vi with discrete fiber of cardinality 2(i, everywhere except on a set of positive 
codimension. To see this, first observe (cf. (|2.2|) ) that the polynomial H has degree 2c? in the 
variable Zd+\. Let "K C V be the subvariety on which dH/dzd+i vanishes. Then on the regular 
set U :— T\ 7r(F), the inverse image of tt contains 2d points and there are distinct, locally defined 
smooth maps yi(x), . . . , j/2ci(x) that are inverted by tt. The fibration 

is the aforementioned fibration with fiber cardinality 2d. 

Next, we apply Cauchy's integral formula with u = — e^+i- Let 5*1 and 5*2 denote the circles in 
of respective radii e"^ and 1 + s, and let Tj := x Sj for j = 1, 2. By ((2?3)) . neither Ti nor T2 
intesects V, so beginning with the integral formula and integrating around Ti, we have 

d+l 

.-r-l , 



1 

27ri 

1 

27ri 



Ti 



d+l 



'F{z) dz 
-''~^F{z)dz 



Ti 



-r-l 



F{z)dz 



T2 



1 

27ri 



d+l 



^F{z)dz. 



Expressing the integral over Tj as an iterated integral over x Sj shows that the quantity in square 
brackets is 



/ / z--iF(z) dzrf+i - / z-'-^F 



(z) dz, 



d+l 



(3.3) 



where Z| denotes (zi, . . . ^Zd). The inner integral is the integral in Zd+i of a bounded continuous 
function of (z|, Zd+i), so it is a bounded function of Z|. We may always write the inner integral as 
a sum of residues. In fact, when Z| G C/ it is the sum of 2c? simple residues, and since Td \ U has 
measure zero, we may rewrite p.3p as 



2'Ki 



■ 2d 
.k=l 



^residue(i^(zt, ■);yk{z^)) 



IZ| 



(3.4) 



On U, we have seen from (|3.1I) that 

RES (Fc?z)(z) = TT* [residue {F{z-^, Zd+i) dz^] (7r(z)) , 
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hence, from the fibration, (13. 4p becomes 

2m I z-'-^RES(i^dz) . 

Because the complement of 7r^^[C/] in Vi has measure zero, we have shown that 

The integral over T2 is 0((1 + s)~^'^)\ because s is arbitrary, sending s — > 00 shows this integral to 
be zero. We have assumed that W is empty, so (13. 5|) becomes the desired conclusion \i.2\ . □ 

The next theorem has the quantum random walk as its main target, however it is valid for a 
general class of rational Laurent series, provided we assume the hypotheses of Lemma [3.21 namely 
torality (j2.3|) and smoothness (W ~ 0). Under these hypotheses, the image of Vi under z ^ (logz)/z 
is a smooth co-dimension-one submanifold M of the flat torus; we let /C(z) denote the curvature 
of M. at the point (logz)/i. Of primary interest is the regime of sub-exponential decay, which is 
governed by critical points on the unit torus. We therefore let K denote the set of directions f for 
which f • X is maximized at x = on the closure of the component of the amoeba complement 
in which we are computing a Laurent series. We also assume ()2.1ip (finiteness of W(f)) for each 
f e K. Observing that z = exp(zx) e W if and only if x is critical for the function r • x on A^, we 
may define r(z) to be the signature of the critical point (logz)/i (the dimension of positive space 
minus dimension of negative space) for the function f • x on . 

Theorem 3.3. Under the above hypotheses, let M he a compact subset of the interior of K such 
that the curvatures /C(z) at all points z € W(f) are nonvanishing for all f G f\f. Then as |r| 00, 
uniformly over f e A/", 

provided that Viog-ff is a positive multiple of r (if it is a negative multiple, the estimate must be 
multiplied by —1). When f ^ K then Or — o(exp(— c|r|)) for some positive constant c, which is 
uniform if r ranges over a compact subcone of the complement o/K. 

Proof: The conclusion in the case where r ^ K follows from Corollary 12.71 In the other case, 
assume r £ M and apply Lemma 13.21 to express Or in the form (|3.2p : 

The chain of integration is a smooth d-dimensional submanifold of the unit torus in M'^+i, so when we 
apply the change of variables z = exp(iy), the chain of integration becomes a smooth submanifold 
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M of the flat torus Tq, hence locally an immersed d-manifold in R''"'"^. We have dz = izdy, so 
F{z)dz/z ~ i'^ F o exp(y) dy and functoriality of RES implies that 

RES — ^ = RES {F o exp dy) . 

After the change of coordinates, therefore, the integral becomes 

d 



I M 

where rj := RES (F o exp dy). By hypothesis, 77 is smooth and compactly supported, so if we apply 
Corollary 12 . 91 and divide by {2ttY we obtain 



1 

27r|r| 



Finally, we evaluate ri[z] in a coordinate system in which the {d + 1)"' coordinate is f . We see 
from ((3?T|l that 

dH/dr{z) 

where dr A dA ~ dz. Because the gradient of H is in the direction f , this boils down to rj = 
G(z)/| Viogi?(z)| at the point z, finishing the proof. □ 



3.2 V contains noncontributing cone points 



In this section, we generalize Theorem [32] to allow VH to vanish at finitely many points of V. The 
key is to ensure that the contribution to the Cauchy integral near these points does not affect the 
asymptotics. This will be a consequence of an assumption about the degrees of vanishing of G and 
H at points of W. We begin with some estimates in the vein of classical harmonic analysis. Suppose 
is a smooth p-form on a smooth cone in ]R''+^; the term "smooth" for cones means smooth except 
at the origin. We say 77 is homogeneous of degree k if in local coordinates it is a finite sum of 
forms A{z) dzi^ A • • • A dzi^ with A homogeneous of degree k — p, that is, A{Xz) ~ A'^~pA(z). A 
smooth p-form 77 on a smooth cone is said to have leading degree a if 



J2 0(|z|"-P+idz,, Adz,J (3.7) 



with 77° hom ogeneous of degree a. The following lemma is a special case of the big-O lem ma 
from BPOSj l . That lemma requires a rather complicated topological construction from ABG70l | ; we 



give a self-contained proof, due to Phil Gressman, for the special case required here. 

Lemma 3.4. Let Vq be a smooth {d— 1)- dimensional manifold in S"^ and let V denote the cone over 
Vq in R'^+^. Let rj he a compactly supported d-form of leading degree a > on V . Then 

e*'-^/ = 0(|r|^"). 
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Proof: Assume without loss of generality that rj is supported on the unit polydisk {z : |z| < 1}, 
where |z| kj P is the usual euclidean norm on C^^^. The union of the interiors of the 

annuli 



< |z| < 2""} 



is the open unit polydisk, minus the origin. Let 0„ ; Bq — > i?„ denote dilation by 2~" and let 
rjn '■— 0'!^f]\Bo be the puUback to Bq from i?„ of the form rj. Let ?f denote the homogeneous part 
of r/, that is, the unique form satisfying (|3.7p . The forms rjn are asymptotically equal to 2~""77° 
in the following sense: for each L, the partial derivatives of 2""77„ up to order L converge to the 
corresponding partial derivatives of 77°, uniformly on Bq. Let x„ be smooth functions, compactly 
supported on the interior of Bq , and with partial derivatives up to any fixed order bounded uniformly 
in n. Then for any N > there is an estimate 



^Xn(z)-(2"'V(z)) = 0(|rr^) 



(3.8) 



uniformly in n. This is a standard result, an argument for which may be found in Ste93l Proposi- 
tion 4 of Section VIII.2], noting that uniform bounds on the partial derivatives of coefficients of XnTJn 
up to a sufficiently high order L suffice to prove Stein's Proposition 4 for the class ?7„, uniformly in 
n. To make the O-notation explicit, we rewrite (|3.8p as 



e''--Xn(z)r/„(z) <<?jv(|r|)2- 



for some functions ^Ar(x) each going to zero as x ^ oo. 



-an |„| — A'' 



(3.9) 



Next, let {ipn : n. > 0} be a partition of unity subordinate to the cover {Bn}- We may choose V'n 
so that < i/",! < 1 and so that the partial derivatives of V'n up to a fixed order L are bounded by 
Cl2" where Cl does not depend on n. We estimate /g e^'^'^^l^nV in two ways. First, using \^pn\ < 1 
and r]{z) = 0{\z\°'~'^ dzi^ ■ ■ ■ dzi^), we obtain 



< C2-"'' sup |z|"-'' < C'2- 



(3.10) 



for some constants C, C independent of n. On the other hand, pulling back by we observe that 
the partial derivatives of O^'^n up to order L are bounded by Cl independently of n. Using p.9p . 
for any iV > we choose L = L{N) appropriately to obtain 



-N 



for all n, N, where gpf are real functions going to zero at infinity. 
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Let no(r) be the least integer such that 2 < l/|r|. Our last estimate implies that for n 

no- j < no, 

-N 



< 2" 



9N{^ If 



2- alio 



2^- 



-N 



2"o / \ 2"" 

Once N > a, the quantity in the square brackets is summable over j > I, giving 



E 

n<no 



Br, 



On the other hand, (|3.10p is summable over n > no, so we have 



E 

n>no 



2- 



The last two estimates, along with |r| = 0(2"°), prove the lemma. 



□ 



Given an algebraic variety V := {H = 0}, let p be an isolated singular point of V. Let H° = H° 
denote the leading homogeneous term of H at p, namely the homogeneous polynomial of some degree 
TO such that H{p + z) = H°{z) + 0{\z\"^~^-^); the degree to will be the least degree of any term in the 
Taylor expansion of H near p. The normal cone to V at p is defined to be the set of all normals 
to the homogeneous variety Vp := {z : H°{p + z) = 0}. We remark that r is in the normal cone to 
V at p if and only if r • z has (a line of) critical points on Vp. 

Theorem 3.5. Let F,G, H,V, Bo,Vi andW be as stated at the beginining of this section. Assume 
torality (j2.3p . Suppose that the singular set W is finite and that for each p G W, the following 
hypotheses are satisfied. 



(i) The residue form r/ has leading degree a > d/2 at p. 

(ii) The cone Vp is projectively smooth and r is not in the normal cone to V at p. 

Then a conclusion similar to that of Theorem \3.3\ holds, namely the sum (j3.6[) over the points Zj ^ W 
where VH \\ r gives the asymptotics of Or up to a correction that is o(|r|^'^/^). 



Proof: By Tou68l . Cor. 2"], condition (ii) implies that the function H{p + z) is bi-analytically 
conjugate to the function H°, that is, locally there is a bi-analytic change of coordinates ^p such 
that H° o v[/p = H{p + z). Now for each p G W, let Up be a neighborhood of p in V sufficiently small 
so that it contains no other p' € W, contains no yj, and so that the bi-analytic map is defined 
on Up. Let Uo be a neighborhood of the complement of the union of the sets Up. Using a partition 
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of unity subordinate to {(7p, C/q}, we replicate the beginning of the proof of Theorem [23] to see that 
it suffices to show 

/ e'^-yRES {F d^) = o{\r\~^/^) . 

Changing coordinates via gives an integral of a smooth, compactly supported form 77 on the cone 
Vp which is homogeneous of order a > d/2. Lemma estimates the integral to be 0(|r|~"), which 
completes the proof. □ 



4 Application to 2-D Quantum Random Walks 

As before, we let F = i<ij<fe where 

r,s,t 

and ai]s^n is the amplitude for finding the particle at location (r, s) at time n in chirality j if is 
started at the origin at time zero in cardinality i. Each entry F^^'^^ has some numerator G(*'^) and 
the same denominator H = det(/ — zMU). In addition, we will denote the image of the gauss map 
of Vi \ W as C/. We note that r E G precisely when 

There is some z in the unit torus for which H{z) = and Viogi? (z) || f . (4-11) 

In fact, we can make a stronger statement as follows. 

Lemma 4.1. G CK. 

Proof of Lemma [4. It Let z satisfy (|4.1ip for some f . Because V is smooth at z, a neighborhood of 
z (or a patch including z) in V is mapped by the coordinatewise Log map to a support patch to Bq 
which is normal to f . This patch lies entirely outside Bq by the convexity of amoeba complements. 
In the limit we see the following. If we take the real version of the complex tangent plane to V G C^^^ 
at z and map by the coordinatewise log map, the result is a support hyperplane to Bq which again, 
lies completely outside Bq (except at Log|z|) by convexity. Now when r G Q, equation (|4.11[) is 
satisfied with z e Vi. Thus Log|z| = and f e K. The desired conclusion follows. □ 

We will apply the results of Section [3] to several one-parameter families of two-dimensional 
QRW's. Each analysis requires us to verify properties of the corresponding family of generating 
functions. 
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4.1 The family S{p) 



We begin by introducing a family S{p) of orthogonal matrices with p G (0, 1): 



Sip) 





VP 


VP 




VT^ 


\ 












VP 


VP 




VT^ 






V2 


V2 










VT^ 


V^ 


VP 








V2 ^ 


v^ 


V2 


V2 




v 






VP 


Vp 


) 


V2 


V2 


V2 


V2 



The matrix S{l/2) is the alternative Hadamard matrix referred to earlier as C/nad- A probability 
profile was shown in figure [21 he re is a picture for another parameter value, namely 1/8. The 
following theorem, conjectured in BraOTj, shows why similarity of the pictures is not a coincidence. 




(a) limit (b) probabilities at time 200 

Figure 5: the S{l/8) walk 



Theorem 4.2. For the quantum random walk with unitary matrix U — S{p), let Q' he a compact 
subset of the interior of Q such that the curvatures /C(z) at all points z G W(f) are nonvanishing 
for all f e Q' . Fix chiralities i,j, let G :— G^^'^\ and let Or := ar,s,n denote the amplitude to be at 
position (r,s) at time n. Then as |r| — s- oo, uniformly over f G Q' , 

«r = (-1)'^ F / + O ( |rr3/2) (4.12) 

27r|r|^^ |V,ogi?(z)| ^^^^ V ' J ^ > 

where 5=1 i/ Viog-ff is a negative multiple of r (so as to change the sign of the estimate) and 
zero otherwise. When f G [— 1, 1]^ \ then for every integer N > there is a C > such that 
Pr (r) < C|r|~^ with C uniform as r ranges over a neighborhood Af of r whose closure is disjoint 
from the closure of Q . 
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Before proving this theorem we interpret its implication for the probabihty profile. The proba- 
bility of finding the particle at (r, s) in the given chiralities at the given time is equal to |arp. We 
only care about Or up to a unit complex multiple, so we don't care whether S is zero or one, but 
we must keep track of phase factors inside the sum because these affect the interference of terms 
from different r e W. In fact, the nearest neighbor QRW has periodicity (because all possible steps 
are odd); the manifestation of this is that W consists of conjugate pairs. When r + s and n have 
opposite parities the summands in the formula for Or cancel. Otherwise the probability amplitude 
jorp will be G(n~^), uniformly over compact regions avoiding critical values in the range of the 
logarithmic Gauss map but blowing up at these values. 

Proof of Theorem 14. 2t As C K by lemma HTTl the result when f e 5' is immediate once we 
have shown that for any S{p), its generating function satisfies the hypotheses of Theorem 13.31 We 
establish this in the lemma below. 

Lemma 4.3. Let H := H^p^ = det {I ~ zM{x,y)S{p)). Then for < p < 1, \JH ^ on T3. 
Consequently, Vi :— Vh H T3 is smooth. 



Theorem l3 . 31 will not be helpful in proving the case when f g [— 1, 1]^ \ ^. To prove this condition 
we present the following lemma, which is a generalization of Ste93l . Proposition 4 of Section VIII.2]. 

Lemma 4.4. Let A4 be a compact d-manifold. Suppose a is smooth and that f is a smooth function 
taking values in R/(2nL), with no critical points in Ai. Then 



/(A) 



(4.13) 



M 



as X 00 through multiples of L, for every N >0. 



We will see below that Vi is a fourfold (unbranched) cover of the two-torus. Any such cover is 
compact. In the calculation of Or, we have /(y) = — f • y and A = |r|. Thus a direction f is not in 
G precisely when /(y) has no critical points in Vi. Uniform exponential decay of amplitudes for r 
bounded outside the image of the gauss map follows. □ 

We now prove the above lemmas in reverse order. 

Proof of Lemma 14.41 : As is compact it admits a finite open cover {C/i}ig/ with subordinate 
partition of unity {(f>i}i,£i. We decompose the integral 
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M 



e 



M 



We will show that for each i e /, Jjj, e^'^^^^'>a{x)<j)i{x)dx is rapidly decreasing (the requirement above 
for /(A)). As the cover Ui is finite, this will give us our result. 

For a given i e /, we let ipi^) ■= ct(2;)0i(x) which is then smooth with compact support. For 
each xq in the support of ip{x), there is a unit vector ^ and a small ball B{xq), centered at xq, such 
that ^- iy f ){x) > c > for some real c uniformly for all x E B{xq). We then decompose the integral 
Su G^^^'"^'^''P{x)dx as a finite sum 

j e'^^^^Hk{x)dx 

k 

where each tpk is smooth and has compact support in one of these balls. It then suffices to prove 
the corresponding estimate for each summand. Now choose a coordinate system xi, . . . ,Xd so that 
xi lies along ^. Then 



, . dXrl 



Now by Ste93l . Proposition 1 of Section VIII. 2] the inner integral is rapidly decreasing, giving us 



our desired conclusion. □ 

For the next two proofs, we clear denominators to obtain the following explicit polynomial: 
H ~ (x^y"^ +y'^ — x^ — l + 2xyz^)z^ — ^xy—y/ljpzixy^ — y — x + z^y — z^x + z'^xy^ + z^x^y — x^y). We 
make the substitution a — \f2/p to facilitate the use of Grobner Bases, which require polynomials as 



inputs. Use the notation Hx for and similarly with y and z. 
Proof of Lemma [4731 

Using the Maple command Basis([H, Hx, Hy, H^], plex(x, y, z, a) we get a Grobner Basis with first 
term za^(a^ — l)(a^ — 2) = 2zp(2p— l)(2p — 2). Thus to show that Siji) results in a variety whose 
intersection with T is smooth for p e (0, 1), we need only consider the case whenp = 1/2. In this case 
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a = 1 and the Grobner Basis for the ideal where {H, Vi7) = is {—z + z^,z^ + 2y~z,—z—z^ + 2x). 
Here Bi vanishes on the unit circle for z = ±1, ±i. However, for z = ±1, B2 vanishes only when 
y = and for z = ±i, B3 vanishes only when x = 0. Thus VTJ does not vanish on T^. □ 

Further analysis of the limit shape for S{p) 

Proposition 4.5. For each pair (x, y), there are four distinct values zi, Z2, Z3, Z4 such that (x, y, Zi) G 
Vi for i e 1,2,3,4. Consequently, the projection {x,y,z) i— > {x,y) is a smooth four-covering of T2 
by Vi. 

Proof: Since H has degree four in z, it has at most four z values in C for each pair (x, y), hence 
at most four z values in Vi- Recall from Proposition 12.11 that all solutions to H{x,y,z) = for a 
given {x,y) in the unit torus have \z\ = 1 as well. Hence, if ever there are fewer than four z values 
for a given (x, y), then there are fewer than four solutions to H{x, y, ■) = and the implicit function 
theorem must fail. Consequently, ^ = 0. This cannot be true, however, by the following argument. 
We have ruled out Hx = Hy = Hz ^ on Vi, so if Hz = 0, then the point {x,y,z) contributes 
toward asymptotics in the direction (r, s, 0) for some (r, s) ^ (0, 0). The particle moves at most one 
step per unit time, so this is impossible. □ 

To facilitate discussions of subsets of the unit torus, we let (a,/3,7) denote the respective argu- 
ments of (x, y, z), that is, x = e*", y = e"^^ , z = e^'^ . We may think of a,/3 and 7 as belonging to the 
flat torus (R/27rZ)3. 

Proposition 4.6. Vi can he decomposed into connected components asVi — AYLBYLCUD, where 
A,B,C and D will be the components containing the 7 values 7:12,1: and 37r/2, respectively. 

Proof: Let x '■= {{x^Viz) : z^ = —1}. We begin by estabhshing that |Vi n x| = 8 with two 
points for each of the fourth roots of —1. Furthemore, — 7r/4 < 7 < 7r/4 on A, nlA < 7 < 37r/4 
on i?, 37r/4 < 7 < 57r/4 on C, and 57r/4 < 7 < Tii/A ox\ D. These observations suffice to prove 
the proposition, because the smooth variety Vi cannot have an intersection with a torus that is 
pinched down to a point; the only possibility is therefore that these values of 7 are extreme values 
on components of Vi . 

Tocheckthe first of these statements, use the identities C0S7 = (z+z^-'^)/2, sin7 = (z — z^^)/(2i), 
as well as the analogous identies for a and /3, to write the equation of V in terms of a, (3 and 7. We 
find that Hix, y, z) = if and only if 

^ L{a, f3,^) -.= 2 sin 7 cos 7 — •\/2p(sin /3 cos 7 -I- cos a sin 7) + cos a sin /? . (4-14) 

Substituting 7 = 7r/4 results in 

1 — (sin /3 + cos a)y/p + cos a sin/3 = . 
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Verifying that sin (3 = ,Jp is not a solution, and dividing by sin [3 ~ ,Jp, wc find that 

1 — sin j3 
cos a = -^-^ — . 

The right-hand side is in [—1, 1] only when sin = ±1. Thus when 7 = 7r/4, the pair (a, /3) is either 
(7r,7r/2) or (0,37r/2). 

To check the remaining statements, we introduce the following set of isometrics for Vi. Define 

0A(a,^,7) := -7) 

/9+ -,7+ -j 

0c(a,^, 7) := (a + TT, /? + TT, 7 + tt) 
(/)D(a,/3,7) := I P + " + Y' Y 

Verifying that 0^1, and (^c (and hence which is equal to (^c ° ^b) are isometrics is a simple 
exercise in trigonometry using equation 14.141 which wc will omit. Each isometry inherits its name 
from the region it proves isometric with A. Using these isometries, we see that 7 is equal to 37r/4, 
57r/4 and 77r/4 exactly twice on Vi. □ 

We remark upon the existence of an additional eight-fold isometry within each connected com- 
ponent: (?!)i(a,^,7) := (a, ^ tt, -7), (l)2{a,/3,j) := {-a,f3,j) and (l>3{a,/3,j) := {a,7T - 
These symmetries manifest themselves in the plots in figures [2] and [5] as follows. The image is 
clearly the superposition of two pieces, one horizontally oriented and one vertically oriented. Each 
of these two is the image of the Gauss map on two of the regions A, B,C, D, and each of these 
four regions maps to the plot in a 2 to 1 manner on the interior, folding over at the bound- 
ary. To verify this, we observe that if po contributes to asymptotics in the direction (r, s) then 
(f>A{po), 4>b{po), 4*0 (po), 4>d{po), 4'iiPo)j 02 (po) and i^siPo) contribute to asymptotics in the directions 
(r, s)(s, r), (r, s), (s, r), (— r, — s), (— r, s) and (r, — s), respectively. Thus while the image of the Gauss 
map is two overlapping leaves, the Gauss map of A and C contribute to one leaf, while the Gauss 
map of B and D contribute to the other. 

We end the analysis with a few observations on the way in which the plots were generated. Our 
procedure was as follows. Solving for sin 7 in (|4.14p . we obtained 



. ^ v2pcos7-cosa -,r\ 

sm 7 = sm /3 -= . (4-15) 

2 cos 7 — V 2p cos a 

Squaring (|4.14p and making the substitution sin'^ 7 = 1 — cos^ 7, we found that 

(1 — cos^ 7) ^2 cos 7 — y/2p cosa^ ~ ~ Z^) cos 7 — cosa^ 

which we used to get the four solutions for 7 in terms of a and (3. We then let a and /? vary over 
a grid embedded in the 2-torus and solved for the four values of 7 to obtain four points in Vi; this 
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Figure 6: The variety Vi for p = 1/2 



is the composition of the first two maps in (jl.ip . Differentiation of iJ(e*", e'^, e'^) = shows that 
the projective direction (r, s,i) corresponding to a point (a,/3, 7) is given by r/t = —d^ / da, s /t = 
—d"f/df3. Imphcit differentiation of (|4.14p then gives four expUcit values for {r/t, s/t) in terms of a 
and (3. This is the composition of the last two maps in (jl.ip , with the parametrization of MP^ by 
{r/t, s/t) corresponding to the choice of a planar rather than a spherical slice. 

4.2 The family A{p) 

We now present a second family of orthogonal matrices A{p) below. In order for the matrices to be 
real, we restrict p to the interval (0, l/\/3). 
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P 
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This family intersects the family S{p) in one case, namely A{l/2) = S{l/2): for any {p,p') G 
(0, 1)^ other than (1/2, 1/2), we have A{p) ^ S{p'). The following theorem follows from Lemma 
along with a new lemma, namely Lemma 14.81 below, analogous to Lemma 14.31 

Theorem 4.7. IfO<p< l/\/3 then Theorem \4-<i\ holds for the unitary matrix A{p) in place of the 
matrix S{p). □ 

Lemma 4.8. Let H := ij'p) = Aet {I ~ zM{x,y)A{p)) . Then for Q < p < l/VS, VH ^ on T3. 
Consequently, Vi := Vh H is smooth. 



Proof of Lemma [4.81 We clear our denominator by setting H := {—xy) * det(/ — MA{p)z), now 
to get 

H = 2{x—l){x+l){y'^ + l)z'^p'^--{—y—x+xy'^+z'^y~x'^y+z'^xy'^ — z'^x+z'^x'^y)zp+{yz'^—x){xz'^+y) . 

As no y^l — p'^ term appears, we can determine a Grobner Basis without making a substitution. 
The Maple command Basis([H, Hx,Hy,Hz],plex(x,y,z,p) delivers a Basis with first term p^z{2p + 
l)(8p2 - 3)(2p2 _ i)(2p - 1). The roots of the first four factors faU outside of our interval (0, l/VS) 
while the root of the last factor corresponds to the matrix S'(l/2) for which we know Vi is smooth 
from the discussion above. □ 

Again we use theorem 13.31 to correctly predict asymptotics for individual directions. We show 
probability profiles for a number of parameter values. 
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Figure 7: The profile for ^4(1/6) shows how QRW approaches degeneracy at the endpoints p — > 0, 1 




Figure 8: p increases from 1/3 to 5/9, switching the direction of the tilt 
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4.3 The family B{p) 

To demonstrate the application of theorem 13.51 we introduce a third family of orthogonal matrices, 
B{p), withpe (0,1). 



We have already seen a walk generated by such a matrix, as Figure [T] depicted the walk generated 
by B{l/2). We note that B{p) is almost identical to S{p) with the one exception being the multi- 
plication of the third row by —1. As was the case with the S{p) walks we can see strong similarities 
between the image of the gauss map and the probability profile for various values of p. 



Figure 9: The image of the Gauss map alongside the probability profile for the B{2/3) walk 

In contrast to the cases of S{p) and A{p), we will not be able to apply Theorem 13.31 because Vi 
is not smooth. 

Theorem 4.9. For the quantum random walk with unitary matrix U = B{p), let Q' be a compact 
subset of the interior of Q such that the curvatures IC{z,) at all points z e W(r) are nonvanishing 
for all f G tj'. Then as |r| — > oo, uniformly over r €i G' , 



When f e [— 1, 1]^ \ then for every integer N > there is a C > such that Pr (r) < C|r| ^ with 
C uniform as r ranges over a neighborhood J\f of r whose closure is disjoint from the closure of Q . 






(4.16) 
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Proof: First, we apply lemma with the lemma being applicable as we will see below that 
Vi := Vh n T3 is a two-fold cover of T2 and thus compact. The conclusion when f £ [—1, 1]^ \ Q 
follows. We get the conclusion in the case where f € ^' by verifying the hypotheses of theorem 13.51 
in the following lemmas. 

Lemma 4.10. Let H := H^p'^ = det {I - zM{x,y)B{p)). Then for < p < 1, the set W = 
{{x, y, z) : (H, VH) = 0} consists only of the four points {x, y, z) — ±(1, 1, \/p/2 ± iy^l — p/2). 

Lemma 4.11. For any < p < 1 we have the following conclusions for each po £ W for the 
generating function associated to the unitary matrix U = B{p). 

(i) The residue form rj has leading degree a > d/2 atpo- 

(ii) The cone Vp^ is projectively smooth and r is not in the normal cone to V at pq. 

Proof of Lemma [4. 101 The proof of Lemma [4. 101 is similar to the corresponding proofs in the two 
previous examples, so we give only a sketch. Computing H from (|2.2p and the subsequent formula 
yields 

H = 2xy{z^ + l)-{x + y + xy'^ + x'^y){z^ + z)^/2p+{Apxy + x^ +x^y'^ + l + y'^)z^ 

= xyz^ ■ [Ap+ (4.17) 
2{z^ + z-^) - {{x + x-^) + [y + y-^)) [z + z'^)^ + {x + x-^){y + y-^)\ , 

Treating p as a parameter and computing a Grobner basis of {H, H.^., Hy, H^} with term order 
plex(x,y,z) one obtains {x^ — x,y — x,z{x'^ — l),z'^ — 2x^/pz + 2x^}. Removing the extraneous 
roots when one of x,y ov z vanishes, what remains is ±(1, 1, z) where z solves — 2y/pz + 2 — 0. □ 

Proof of Lemma [4. Ill Condition (z) follows from the fact that for eachpo S W, the denominator 
G'-P-' (x, y, z) vanishes as well as the numerator H^^^ which only vanishes to order 1. To prove (zi), we 
compute the local geometry of {H — 0} near the four points found in the previous lemma. We will 
do this for the points with positive (x, y) ~ (1,1); the case {x,y) — (—1,-1) is similar. Substituting 
X = 1 + u,y — 1 + v,z = zq + w into H and then reducing modulo Zg — 2y^2;o + 2, we find that 
the leading homogeneous term in the variables {w, v, w\ is 4[y/p(1 — p){u^ + ) - {2-p)w\ For 
< p < I, this is the cone over a nondegenerate ellipse and therefore smooth. The dual cone is the 
set of (r, s, t) with + s"^ ~ i'^-ii)^ ^'^ ' '^^^ minimum value of on [0, 1] is greater than 4, 

while the vectors (r, s, t) inside the image of the Gauss map all have + < 4t^, whence r is never 
in the normal cone to V at po- D 

Beginning with (I4.18p . we see that (x, y, z) e Vi 

2 cos^ 7 — (cos a + cos/3) -\/2pcos7 + cos a cos /? + p — 1 = . (4.18) 
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Thus for given a and (3, the four values of 7 are given explicitly by 
7 = ± arccos 



(cos a + cos (3) -J2p ± \ / 2p (cos a + cos f3) — 8 cos a cos (3 



(4.19) 



We then differentiate 14.181 with respect to a and (3 to obtain the partial derivatives 

(?7 sin a cos a — cos 7 



and 



9a sin 7 cos a + cos /? — 4 cos 7 
97 sin /3 cos a — cos 7 



d/3 sin 7 cos a cos /3 — 4 cos 7 
Remark. The fact that we can solve explicitly for Z with this family allows us to more clearly depict 
the connection between curvature and asymptotics. Using Proposition 12 . 31 and (|4.19p . we let Maple 
evaluate V as well as 



dpda 



We then plot JC against — and — ^ as (a, (3) varies over the two-dimensional torus. 




Figure 10: A graph of curvature versus direction for the B{\/2) walk 
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In the above picture we see the expected cross within a diamond region where curvature is low, 
though the view is obstructed by regions of higher curvature. 

To remedy this problem we restrict our view of the /C axis to focus on the smallest values of K, 
which in turn contribute to the largest probabilities. The resulting picture thus predicts the regions 
that will appear darkest in the probability profile. 




^ I I I I I I I I I I I I I I I I I I I I I I I r 
-0.6 -0.35 -0.1 0.15 0.4 



Figure 11: A graph of the aread of lowest curvature and hence highest probabilites for the B(l/2) 
walk 
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